АППРОКСИМАЦИЯ СИГНАЛОВ ПРИ ПОСТРОЕНИИ АВТОМАТИЧЕСКИХ СИСТЕМ 
ОБНАРУЖЕНИЯ УТЕЧЕК НА НЕФТЕПРОВОДАХ 

В. А. Якупов, А.М. Кориков* 

Томский политехнический университет 

Томский государственный университет систем управления и радиоэлектроники 
Е-таіІ: ]ѵа@еІе5у.ги, когікоѵ@а$и.Іи5иг.ги 

Рассмотрены проблемы построения систем автоматического обнаружения утечек, предложены варианты решения задачи до- 
ставки обрабатываемых данных в центры их сбора и обработки. Показана целесообразность использования аппроксимации 
данных для сжатия и фильтрации текущих измерений. 


Введение 

Развитие транспортной сети нефте- и нефте- 
продуктопроводов (НПС) сопровождается непре- 
рывным повышением уровня их технической осна- 
щенности, широким внедрением современных 
средств компьютеризации и автоматизации про- 
цессов перекачки. Трубопроводная транспорти- 
ровка неразрывно связана с выполнением ком- 
плекса мероприятий по охране окружающей среды 
на принципиально новых научно-технических ос- 
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новах проектирования, строительства и эксплуата- 
ции магистральных трубопроводов. Несмотря на 
очевидные преимущества трубопроводного транс- 
порта при транспортировании нефти и нефтепро- 
дуктов по магистральным трубопроводам, имеет 
место загрязнение атмосферы, водоемов и почвы. 
В связи с этим серьезное внимание уделяется со- 
кращению потерь на магистральных трубопрово- 
дах. Как правило, потери связаны с нарушением 
правил эксплуатации, повреждением трубопрово- 
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дов от коррозии, несвоевременным ремонтом, сти- 
хийными бедствиями и участившимися случаями 
несанкционированных врезок в трубопровод с це- 
лью хищения нефти и нефтепродуктов [1]. В связи 
с этим особое значение приобретает внедрение ав- 
томатических систем обнаружения утечек в систе- 
мы диспетчерского контроля трубопроводной 
транспортировки. 

О проблемах построения автоматических 

систем обнаружения утечек 

В настоящее время для обнаружения несанкци- 
онированных врезок и утечек на НПС разработано 
большое количество методов [2], основанных на 
различных физических законах и явлениях. 

Среди многообразия предлагаемых методов, 
наибольший интерес представляют методы обна- 
ружения утечек, пригодные для создания автома- 
тических систем обнаружения утечек. Например, 
лазерный газоаналитический, ультразвуковой зон- 
довый, магнитные, тепловые, радиоактивные, ви- 
зуальные методы неудобны для построения авто- 
матических систем обнаружения утечек, т. к. явля- 
ются в основном методами активного контроля с 
привлечением дорогостоящих материалов и 
устройств. Кроме того, у многих из них имеются 
погодные ограничения, ограничения по глубине 
залегания трубопровода и ограничения по приме- 
нимости в разные времена года. 

Основными требованиями, предъявляемыми к 
системам автоматического контроля НПС, явля- 
ются: высокая чувствительность; точность опреде- 
ления места утечки; безопасность в эксплуатации; 
обеспечение контроля трубопроводов большой 
протяженности; высокая степень надежности, до- 
стоверности и автоматизации; экономичность; го- 
товность к работе при любых климатических и по- 
годных условиях; оперативность предоставления 
результатов анализа данных. 

Одним из наиболее пригодных методов для соз- 
дания автоматической системы непрерывного кон- 
троля герметичности трубопровода является метод 
отрицательных волн давления. В основе метода ле- 
жит регистрация фронта волны изменения давле- 
ния, которая возникает в месте разгерметизации 
трубопровода. По разнице времени прихода фрон- 
та волны на концы контролируемого участка рас- 
считывается место возникновения утечки. 

Данный метод предъявляет высокие требования 
к измерительным системам и средствам связи. В 
ряде случаев в силу объективных причин суще- 
ствующими средствами телемеханики удовлетво- 
рение предъявляемых требований не всегда воз- 
можно. Поэтому не может быть решена не только 
задача обнаружения утечек малой интенсивности 
(до 3 % от текущего расхода) на трубопроводе, но и 
задача адекватного контроля давления с точки зре- 
ния определения критичных значений технологи- 
ческого процесса. Эта проблема характерна для 


всех методов автоматического непрерывного кон- 
троля герметичности трубопровода, которые для 
получения более точного ответа используют боль- 
шое количество информации. 

Наиболее узким местом в этой проблеме явля- 
ется строгая иерархическая система существующей 
телемеханики, при которой данные собираются и 
обрабатываются единым центром. Это приводит к 
ограничению объемов информации, которая мо- 
жет быть передана к центру, с ростом числа узлов и 
количества данных для обработки особенно при 
низкой пропускной способности линии связи. 

Одно из решений устранения данного недостат- 
ка, не приводящее к большим капитальным затра- 
там, является частичная обработка получаемых 
массивов данных на месте, что достаточно эффек- 
тивно решается сегодня путем использования ми- 
кропроцессорных средств, в частности процессо- 
ров обработки сигналов и управляющих контрол- 
леров. При этом диспетчер получает данные, толь- 
ко из выбранного частотного диапазона и времен- 
ного интервала работы, в то время как весь объем 
информации обрабатывается предварительно на 
месте самим контроллером. Остальные данные мо- 
гут накапливаться в энергонезависимой памяти и 
при необходимости быть считаны дистанционно 
или обслуживающим персоналом. 

Вопрос обнаружения утечки на фоне техноло- 
гических помех может быть решен путем объедине- 
ния датчиков давления в систему, при которой ди- 
намика давления будет сопоставляться на различ- 
ных контролируемых участках. Тогда при соответ- 
ствующей синхронизации имеется возможность 
определения местоположения источника всплеска 
давления и, соответственно, его отнесения к тому 
или иному типу. 

Сложность реализации такого подхода заклю- 
чается в том, что технические характеристики вы- 
числительной техники контроллерного уровня су- 
щественно ниже характеристик ЭВМ диспетчер- 
ского уровня, что не позволяет применять эффек- 
тивные и в тоже время затратные алгоритмы обра- 
ботки данных. Поэтому в качестве второго реше- 
ния поставленной проблемы предлагается сокра- 
тить объем передаваемой информации в централь- 
ный пункт сбора и обработки. 

В качестве пути решения задачи сокращения 
объема передаваемых данных предлагается ап- 
проксимация сигнала, представленного массивом 
данных, на контроллерном уровне и передачи ко- 
эффициентов аппроксимации в центр сбора и об- 
работки данных. Такой подход позволит использо- 
вать простые и малозатратные алгоритмы аппрок- 
симации на местах и эффективные алгоритмы об- 
наружения утечек на диспетчерском уровне. 

Аппроксимация текущих измерений 

Основные трудности при решении задачи ап- 
проксимации состоят в следующем: выбор базиса, 
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позволяющего представить данные требуемого ча- 
стотного диапазона; малое количество коэффици- 
ентов для представления данных с допустимой по- 
грешностью; реализация алгоритма, работающего 
в режиме реального времени. 

Для систем, функционирующих в реальном 
времени, важна простота реализации на ЭВМ, по- 
этому особый интерес представляет полиномиаль- 
ная аппроксимация для кусочно-постоянных пара- 
метров. В данной работе предлагается использова- 
ние кубических сплайнов дефекта равного двум в 
качестве базиса алгоритма аппроксимации. В поль- 
зу такого решения говорят следующие факты: для 
кубического сплайна с требованием сопряжения 
по значению и по первой производной необходимо 
рассчитать только два коэффициента. Существует 
несколько реализаций алгоритма, имеющих в ос- 
нове итерационный подход и способных функцио- 
нировать в режиме реального времени. К таким ал- 
горитмам можно отнести, например, построение 
рекуррентного сглаживающего сплайна на основе 
метода наименьших квадратов [3] или применение 
алгоритма Гревилля решения систем линейных ал- 
гебраических уравнений при постановке задачи в 
терминах псевдообратных матриц [4]. 

Рассмотрим подробнее задачу нахождения ко- 
эффициентов кубического сплайна. Алгебраиче- 
ский полином представляется в следующем виде: 

Р(1) = с/ +с п _/~' + ... + с/+с 0 , (1) 

поэтому определение коэффициентов с 0 ,с 1ѵ ..с„ сво- 
дится к решению системы линейных алгебраиче- 
ских уравнений (СЛАУ) вида 

С„і" + С п _ х 1"~' + ...+ с 1 1 1 +с 0 ау і , 1 = 1 , т. ( 2 ) 


Для удобства последующего изложения введем 
в рассмотрение матрицу А и вектор-столбцы сиу, 
определив их равенствами: 



(и і 2 

і п л 




м 

А = 

1 с с 

і п 

, . і 2 

4. 

,с = 

с. 

4 

,У = 



У*»*» • 

і п , 

•• 1 т ) 


V е » У 


У т ) 


Представим систему (2) в векторно-матричной 
форме 
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В общем случае для задачи аппроксимации 
представленная система (3) будет несовместной. 
Несовместность системы уравнений (3) обуславли- 
вает необходимость обобщения классического по- 
нятия «решение СЛАУ». Подход к введению по- 
добного обобщения заключается в том, что задача 
решения СЛАУ сводится к решению экстремаль- 
ной задачи: минимизация выбранной количе- 
ственной меры погрешности аппроксимации. По- 
скольку решение экстремальных задач существует 
независимо от того, является ли система (3) совме- 
стной или нет, то реализация данного подхода по- 


зволяет полностью решить проблему существова- 
ния решения СЛАУ. При таком подходе, учитывая 
свойства псевдорешения СЛАУ [5], целесообразно 
использовать ее псевдорешение, определяемое ра- 
венством вида 

4 

с = А* у. 


где А + - матрица, псевдообратная к матрице А. 

Среди известных алгоритмов вычисления ма- 
трицы А* в данной работе предпочтение отдается 
алгоритму Гревилля. Выбор этого алгоритма целе- 
сообразен из-за его рекуррентной структуры, кото- 
рая позволяет регулировать число вычислений, со- 
гласно выбранной точности аппроксимации. Для 
удобства описания выбранного алгоритма будем 
рассматривать матрицу А как совокупность вектор- 
столбцов и введем в рассмотрение последователь- 
ность матриц А р определяемых равенствами: 
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где й у — у-й столбец матрицы А, /— І ,п. 

Из этих равенств видно, что каждая из матриц А і 
составлена из первых у столбцов матрицы А и связа- 
на с пре дш ествующей ей матрицей соотношением 
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где вертикальная черта отделяет последний стол- 
бец матрицы. 

Используя приведенную выше последователь- 
ность матриц Ар алгоритм Гревилля можно пред- 
ставить в виде следующей последовательности опе- 
раций: 
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где Т- операция транспонирования; (.,.)- опе- 
рация скалярного умножения векторов; 0,„ - нуле- 
вой да-мерный вектор-столбец; Ь і - вектор-строка 
размерности /я; А* - псевдообратная к А і матрица. 

Сущность алгоритма решения задачи аппрокси- 
мации таблично заданной функции и реализующих 
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его операций изложены применительно к случаю, 
когда в качестве аппроксимирующей функции ис- 
пользуется алгебраический полином вида (1). 

Следует учитывать, что механизм получения 
псевдорешений оправдан в случае неравномерной 
сетки. В случае равномерной сетки процедура по- 
лучения псевдообратной матрицы существенно 
упрощается. Рассмотрим это на примере для куби- 
ческого сплайна дефекта 2. Решение ищется в виде 

4. Г 
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- вектор измерении. 


Так как временная сетка равномерна, т. е. 
|г,_і — /,|=АГ, то с учетом локализации временного ин- 
тервала (С = 0) матрица А принимает следующий вид: 
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Из условий сопряжения сплайна и его первой 
производной (с 0 и С[ известны) получаем систему, 
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то есть вычисление псевдообратной матрицы сво- 
дится к определению Известно что 

X ' 4 = 4 т л ( 2 и - !)(« - !)( 3 « 2 - 3 » - 1 ); 

І=1 ^ 

|>' 5 =і» 2 (2« 2 -2п-1)(п-1) 2 ; 

; 6 = -^гп(2п - \)(п - 1)(3я 4 - 6 п + 3 п + 1). 

і=і ^2 

Нетрудно видеть, что число вычислений при 
определении / не изменяется при росте числа изме- 
рений, участвующих в аппроксимации, что нема- 
ловажно для систем, функционирующих в реаль- 
ном масштабе времени. 

Определение размера массива измерений, под- 
лежащих аппроксимации, определение точки со- 
пряжения сплайнов обусловлено требованиями за- 
дачи и выбирается разработчиком исходя из этих 
требований. Рассмотрим эту задачу на примере 
применение кубического сплайна для аппроксима- 
ции текущих измерений давления для последую- 
щей передачи информации по каналам связи и об- 
работки. 

Для обнаружения утечек малой интенсивности 
и определения координаты утечки с погрешностью 
не более 1 % от длины контролируемого участка 
необходимо измерение давления с дискретностью 
порядка 100... 120 Гц, разрядность цифро-аналого- 
вого преобразователя не менее 24 бит. Организация 
файлового обмена (не протоколируемые данные 
телемеханики) между контроллерным уровнем и 
центром сбора и обработки данных позволит со- 
кратить трафик за счет уменьшения избыточности 
универсальных протоколов передачи данных. Для 
метода отрицательных волн давления желательно 
наличие двух датчиков на каждом контрольном 
пункте для определения направления волнового 
возмущения, что увеличит объем передаваемых 
данных с одного контрольного пункта в два раза. 
Поэтому, если выбрать файловый интервал 180 с и 
частоту дискретизации 100 Гц, получим, что объем 
передаваемых данных (без учета сопроводительной 
информации) составит около 105 Мб. 

Выберем число текущих аппроксимируемых 
данных равное 200 и потребуем выбор точки сопря- 
жения в середине сплайна. При аппроксимации 
данных кубическими сплайнами, с условием со- 
пряжения по значению и первой производной, при 
расчете 4 байта на один коэффициент, получим 
5,6 Мб. Если учесть, что первые два коэффициента 
сплайна при меньших степенях рассчитываются по 
коэффициентам предыдущего сплайна, то в файле 
достаточно передать только один раз четыре коэф- 
фициента, а затем передавать только «старшие» ко- 
эффициенты, поэтому объем передаваемых данных 
сократится до 2,8 Мб. 

Таким образом, варьируя количеством аппрок- 
симируемых данных, можно получить требуемый 
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Рисунок. Частотные характеристики сплайн-преобразования 


коэффициент сжатия и передавать на обработку 
максимально возможное количество данных с уче- 
том погрешности, вносимой алгоритмом аппрок- 
симации. 

Для представления частотных характеристик 
сплайнового преобразования с указанными пара- 
метрами воспользуемся средствами среды МАТ- 
ІАВ 6.5 [6]. Частотные характеристики, предста- 
вленные на рисунке, свидетельствуют о том, что 
аппроксимированный сигнал содержит информа- 
цию о волновых процессах с частотой 1 Гц и ниже, 
что соответствует волновым процессам, возника- 
ющим при разгерметизации трубопровода. 

Очевидно, что с изменением количества точек 
аппроксимации и выбора точки сопряжения изме- 
нятся частотные характеристики сплайн-преобра- 
зования. 


Заключение 

Выбор функциональной структуры автоматиче- 
ской системы обнаружения утечек является важным 
моментом при решении задач начального этапа раз- 
работки, который предопределит дальнейшее разви- 
тие системы. Предлагаемый подход предваритель- 
ной обработки данных на контроллерном уровне по- 
зволяет построить систему, в которой основная рабо- 
та по их анализу выполняется на диспетчерском 
уровне. Такое разделение функциональной части си- 
стемы удобно для применения эффективных и в то- 
же время затратных алгоритмов обработки и анализа 
данных. В пользу предлагаемого подхода говорит 
также тот факт, что гибкость систем верхнего уровня 
позволяет в короткие сроки проводить модифика- 
цию и масштабирование, что немаловажно вслед- 
ствие затратности процесса интегрирования систем. 
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